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^ . Abstract 
> 

\jD I Originally designed for CP violation studies in the B meson system, 

^ ' the B-Factories recently showed an exciting capability for improving our 

Ti^ij- ■ experimental knowledge in the field of hadron spectroscopy. We review 

. results on bottomonium spectroscopy from the Belle experiment at the 

O \ KEK-B e"'"e~ collider and present exciting new results from the unique 

large data set taken at the T(55') resonance. 
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1 Introduction 



The Belle Collaboration has collected a large sample of e~^e~ collisions at the en- 
ergy of the T{5S) resonance, which lies above the threshold for production of Bg 
meson pairs, primarily with a purpose of studying decays of Bg. There have been 
a number of unexpected results on the non-BsBg decays of the T{5S). In partic- 
ular, anomalously large rates for dipion transitions to lower bottomonium states 
T{5S) — (T(15'), T(25'), T(35'))7r"''7r~ have been observed [Tj. If these signals are at- 
tributed entirely to the T{5S) decays, the measured partial decay widths r[T(5S') — )• 
T(?T,S')7r"^7r~] ~ 0.5 MeV are about two orders of magnitude larger than typical widths 
for dipion transitions among T{nS) states with n < 4. 

Recently the CLEO-c Collaboration observed the process e"*"e~ — )■ hc{lP)'n'~^7T~ at 
a rate comparable to the process e+e~ — )■ J/ipn^n^ at ^/s = 4170 MeV and found an 
indication of an even higher transition rate at the y(4260) energy [2]. This implies 
that the hb{mP) production might be enhanced in the region of the Yf, and motivates 
a search for the hf,{mP) in the T{5S) data. 

We use the full T(5S) data sample with the integrated luminosity of 121.4 fb~^ 
collected near the peak of the T(55') resonance with the Belle detector [3] at the 
KEKB asymmetric-energy e~^e~ collider [1]. 

2 Observation of the hb{lP) and hb{2P) 

We observe the /if,(lP) and /if,(2P) in the missing mass spectrum of tt+tt^ pairs. The 
7r+7r~ missing mass is defined as MM(n^n^) = ^y{Ec.m. — E*+^_y — where 
Ecm. is the center-of-mass (cm.) energy, E*+^_ and are the vr+vr" energy and 

momentum measured in cm. frame. The details of the analysis can be found in [5]. 
The MM(7r"'"7r~) distribution for the selected 7t~^7i~ pairs is shown in Fig. [It^a). In 
this figure only the T(5S') — )■ T(lS')7r+7r^ and T(5S') — )■ T(2S')7r+7r~ transitions are 
discernible. 

To fit the MM{tt^7t~) spectrum, we separate it into three adjacent regions with 
boundaries at MM(7r+7r-) = 9.3 GeV/c^, 9.8 GeV/c^, 10.1 GeV/c^ and 10.45 GeV/c^. 
We fit every region separately to better control the complicated shape of the combi- 
natorial background, which is described by a Chebyshev polynomial of 6- 7th order. 
In region 3 we subtract the Kg contribution bin-by-bin, while in other regions its 
shape is smooth and is absorbed into combinatorial background. The signal peaks 
are described by Gaussians with paramenters obtained from exclusive decays of the 
T[nS) to ii~ ■ The MM{tt^tt~) spectrum with the combinatorial background and 

contributions subtracted, and the signal function resulting from the fit overlaid, 
are shown in Fig. [U^b). The significance of the hb{lP) and hh{2P) signals which 
includes the systematic uncertainty is 5.5a and 11.2cr, respectively. 
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Figure 1: (a)The MM{7t^7t~) distribution for the selected 7r~^7r~ pairs. Vertical lines 
indicate the locations of the T(IS'), hb{lP), T{2S), hb{2P) and T{3S) signals, (b) 
The MM{'ir~^7c~) spectrum with the combinatorial background and Kg contributions 
subtracted (dots with error bars) and signal component of the fit function (solid 
histogram). The vertical dashed lines indicate the boundaries of the fit regions. 



This is the first observation of the hh{lP) and hi,{2P) spin-singlet bottomo- 
nium states in the reaction e^e" — )■ hbimP)^^^^ at the T(55') energy. We mea- 
sure the masses and the cross sections relative to the e~^e~ — )■ T(25')7r^7r~ cross- 
section: M = 9898.25 ±1.06t}:E]^MeV/c2, R = 0.407 ± 0.079l[j:[}^^ for the h{lP) and 
M = 10259.76 ± OMtlilMeV/c^, R = 0.78 ± OM^ofo for the hb{2P). The masses 
do not differ significantly from the center-of-gravity of the corresponding Xbj states. 
For the hyperfine splitting we find AMhf = 1-62 ± 1.52MeV/c^ for the hb{lP) and 
OAStlilMeV/c^ for the /ib(2P). 

The values of R comparable with unity indicate that the hb{lP) and hb{2P) are 
produced via an exotic process that violates the suppression of heavy quark spin- 
flip. For further study we investigate resonant substructure of these decays Be- 
cause of high background Dalitz plot analysis is impossible with current statistics, 
therefore we study the one- dimensional distributions in M{hb{mP)7r). We define the 
M[hb{mP)TT^) as a missing mass of the opposite-sign pion, MM{tt~). We measure 
the yield of signal decays as a function of the MM{t[^) by fitting the MM(7r"'"7r^) 
spectra in the bins of MM[-n^). We combine the MM{t[^ti^) spectra for the corre- 
sponding MMiiT^) and MM{tt~) bins and we use half of the phase space to avoid 
double counting. 

Results of the fits for the hb{lP) yield as a function of MM{'k) are shown in 
Fig. [21 The hb{lP) yield exhibits a clear two-peak structure without any significant 
non-resonant contribution. In the following we refer to these structures as ^^(lOGlO) 
and Z;,( 10650), respectively. 

We perform a fit to the MM{tt) distributions. We assume that spin-parity 
for both Zf,(10610) and ^^(10650) is = 1"*", therefore in the fit function we use a 
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Figure 2: Left: the yield of the hb{lP) as a function of MM{n) (points with error 
bars) and results of the fit (histogram). Right: the yield of the hh{2P) as a function 
of MM{7t) (points with error bars) and results of the fit (histogram). 



coherent sum of two P-wave Breit-Wigner amplitudes; we add also a non-resonant 
contribution. 

f = A \BW{s, Ml, Ti) + ae''''BW{s, M2, T2) + be''f'\^ ^. (1) 

Here y/s = MM{tt); the variables A, M^, Tk {k = 1,2), a, (j), b and -ip are floating 
in the fit; ^ is a phase-space factor, p [q) is the momentum of the pion originating 
from the T(55') (Zh) decay measured in the rest frame of the corresponding mother 
particle. The results of the fit are shown in Fig. |2] and are summarized in Table [H 
The non-resonant amplitude is found to be consistent with zero. We find that the 
hypothesis of two resonances is favored over the hypothesis of a single resonance (no 
resonances) at the 7.4 a (17.9 a) level. The parameters of the ^^(10610) and Zfe(10650) 
obtained in the fit of hb{lP) and hb{2P) are consistent with each other. 

3 Analysis of T(55) T(15, 25, 35)7r+7r- 

To select T{5S) — )■ T(r2S')7r+7r^ candidate events we require the presence of a pair of 
muon candidates with an invariant mass in the range of 8.0 GeV/c^ < M(/i+/i~) < 
11.0 GeV/c^ and two pion candidates of opposite charge. These tracks are required 
to be consistent with coming from the interaction point. We also require that none 
of the four tracks be positively identified as an electron. No additional requirements 
are applied at this stage. 

Candidate T{5S) — > T(nS')7r+7r^ events are identified by the invariant mass of 
the n~^fJ^~ combination and the missing mass MM{n~^n~) associated with the vr+vr^ 
system calculated as MM(7r~^7i~) = ■>/ (-Ec.m. — ^I+tt- )^ ^ P*n+TT- ' where -Ec.m. is the 
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Figure 3: Dalitz plots for T(nS')7r+7r~ events in the (a) T(IS'); (b) T{2S); (c) T(3S') 
signal regions. Dalitz plot regions to the right of the vertical lines are included in the 
amplitude analysis. 

center-of-mass (cm.) energy and E*+^^ and are the energy and momentum of 

the 7r"^7r~ system measured in the cm. frame. 

The amplitude analyses of the three-body T{5S) — )■ T(nS')7r+7r~ decays that are 
reported here are performed by means of unbinned maximum likelihood fits to two- 
dimensional Dalitz distributions. 

Before fitting the Dalitz plot for events in the signal region, we determine the 
distribution of background events over the Dalitz plot using events in the T{nS) 
mass sidebands that are refitted to the nominal mass of the corresponding T{nS) 
state to match the phase space boundaries. 

In the sideband Dalitz distributions one can see a strong concentration of back- 
ground events in the very low vr^vr^ invariant mass region; these are due to photon 
conversion on the innermost parts of the Belle detector. Because of their low energy, 
these conversion electrons are poorly identified and pass the electron veto require- 
ment. We exclude this high background region by applying the requirements on the 
7r+7r~ invariant mass. For the remainder of the Dalitz plot the distribution of back- 
ground events is assumed to be uniform. The variation of reconstruction efficiency 
across the Dalitz plot is determined from MC simulation. The fraction of signal events 
in the signal region for each of the three T(n5')7r+7r~ final states is determined from 
a fit to the corresponding MM {7r~^7r~) spectrum using a Crystal Ball function [7j for 
the T signal and a linear function for the combinatorial background component. 

Figure [3] shows Dalitz plots of the events in the signal regions for the three decay 
channels under study. In all cases, two horizontal bands are evident in the T[nS)7r 
system near 10.61 GeV/c^ (~ 112.6 GeVVc^) and 10.65 GeV/c^ (~ 113.3 GeVVc^). 

We use the following parameterization for the T{5S) — T(n5')7r+7r~ three-body 
decay amphtude: 

M(si, S2) = Ai{si, S2) + ^2(51, S2) + + Af^ + Anr, (2) 



4 



where si = M^{T{nS)TT+), S2 = M'^(T(nS)'K ). Here we assume that the dominant 
contributions come from the amphtudes that conserve the orientation of the spin of the 
heavy quarkonium state and, thus, both pions in the cascade decay T{5S) — > ZbH — t- 
T(?7.S')7r"'"7r~ are emitted in an S-wave with respect to the heavy quarkonium system. 
As will be shown later, angular analyses support this assumption. Consequently, 
we parameterize the observed Zb(10610) and ^^(10650) peaks with an S'-wave Breit- 
Wigner function without s dependence of the resonance width F. To account for the 
possibility of T(55') decay to both Z^n^ and Z^n^, the amplitudes Ai and A2 are 
symmetrized with respect to 7r+ and tt~ transposition. Taking into account isospin 
symmetry, the resulting amplitude is written as 

A,{si,S2) = a,e''''{BW{si,M,,T,) + BW{s2,Mk,Tk)), (3) 

where the masses Mk and the widths {k = 1,2) are free parameters of the fit. Due 
to the very limited phase space available in the T{5S) — >■ T(3S')7r+7r~ decay, there is 
a significant overlap between the two processes T{5S) — )■ and T{5S) — )■ Z^7r+. 

We also include amplitudes Aj^ and Af^ to account for possible contributions in the 
TT+TT" channel from /o(980) scalar and /2(1270) tensor states, respectively. Inclusion 
of the /o(980) state is necessary in order to describe the prominent structure in the 
M{7r~^7i^) spectrum for the T(lS')7r+7r^ final state around M(7r^7r~) = l.OGeV/c^ 
(see Fig. H]). We also find that the addition of the /2(1270) gives a better description 
of the data at M(7r+7r^) > l.OGeV/c^ and drastically improves the fit likelihood 
values. We use a Breit-Wigner function to parameterize the /2(1270) and a coupled- 
channel Breit-Wigner (Platte) function [8J for the /o(980). The mass and the width 
of the /2(1270) state are fixed at their world average values [9]; the mass and the 
coupling constants of the /o(980) state are fixed at values defined from the analysis 
of B+ K+n+n--. M(/o(980)) = 950 MeV/c^, g^^ = 0.23, qkk = 0.73 [lOj. 

Following suggestions given in Refs.[llj and references therein, the non-resonant 
amplitude A^r has been parameterized as 

^NR = ar-e*'^i" + 4'--e'^?'^-S3, (4) 

where S3 = M^(7r+7r^) (53 is not an independent variable and can be expressed via 
Si and S2 but we use it here for clarity), a^, (12'^, S'^^ and 5" are free parameters of 
the fit (with an exception of the T(35')/T+7r^ channel as described below). 
The logarithmic likelihood function C is then constructed as 

C = -2J2 log(/sig5(si, S2) + (1 - /.ig)5(si, S2)), (5) 

where S{si, S2) = \M{si, 52)^ folded with the detector resolution function (5.6 MeV/ 
for M(T(nS')7r^); the M(7r+7r~) resolution is better and is not taken into account since 
no narrow resonances are observed in the tt+tt" system), B{si,S2) = 1 and /sig is a 
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Figure 4: Comparison of fit results (open liistogram) witfi experimental data (points 
with error bars) for events in the T(IS'), T{2S) and T(35') signal regions. The hatched 
histogram shows the background component. 

fraction of signal events in the data sample. Both S{si, S2) and B{si, S2) are corrected 
for reconstruction efficiency. 

In the fit to the T(15')7r^7r~ sample, the amplitudes and phases of all of the com- 
ponents are allowed to float. However, in the cases of T(2S')7r+7r^ and T(3S')7r+7r^ 
the available phase space is significantly smaller and contributions from the /o(980) 
and /2(1270) are not well defined. Thus, in the fit to the T(25')7r+7r~ and T(35')7r+7r~ 
signal samples, we fix the amplitudes and relative phases of these components to the 
values measured in the fit to the T(lS')7r+7r~ sample. Moreover, in the fit to the 
T(3S')7r+7r^ sample, we also fix the Og'' and ^2"^ parameters of the A^,. amplitude. 
Possible effects of these assumptions are considered while determining the model- 
dependent uncertainty. Results of the fits to T{5S) — )■ T(nS')7r"'"7r~ signal events 
are shown in Fig. HJ where one-dimensional projections of the data and fits are com- 
pared. To combine and events we plot T(nS')7r mass distributions in terms 
of M(T(n5')7r)mm and M{T{nS)n)^ax] fits are performed in terms of M(T{nS)'K^) 
and M{T{nS)7c^). Results of the fits are summarized in Table [H We try various al- 
ternative models to parameterize the decay amplitude as described in the systematic 
uncertainty section. The combined statistical significance of the two peaks exceeds 
10 sigma for all tested models and for all T(n5')7r"''7r~ channels. 

4 Discussion and Conclusions 

In conclusion, we have observed two charged bottomonium-like resonances, the Zfe(10610) 
and Zfe(10650), with signals in five different decay channels, T(?T,S')7r^ (n = 1,2,3) 
and hb{mP)7r^ {m = 1,2). Parameters of the resonances as measured in different 
channels are summarized in Table [TJ All channels yield consistent results as can be 
seen in Fig. [51 A simple weighted averages over all five channels give M[Zfe(10610)] = 
10608.4 ± 2.0 MeV/c2, r[Zfe(10610)] = 15.6 ±2.5 MeV and M[Zf,(10650)] = 10653.2 ± 
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Table 1: Comparison of results on Zfc(10610) and Z;,(10650) parameters obtained from 
T{5S) T(n5)7r+7r- (n = 1,2,3) and T(55) ^ hb{mP)TT+Tr- (m = 1,2) analyses. 
Quoted values are in MeV/c^ for masses, in MeV for widths and in degrees for the 
relative phase. Relative amplitude is defined as aZb(io650)/'^Zi,io6io- 



Final state 


T(15)7r+7r- 


T(2S)7r+7r- 


T{35)7r+7r- 


hb(lP)n+TT- 


hti2P)Tr+n- 


A/(Z6(10610)) 
r(Zi,(10610)) 
M(Zb(10650)) 
r(Zj, (10650)) 
Rel. amplitude 
Rel. phase, 


10609 ± 3 ± 2 
22.9 ± 7.3 ± 2 
10660 ± 6 ± 2 
12± 10± 3 
0.59±0.19lg°^ 
53±6lt^o' 


10616 ± 2ti 
21.1 ±4^3 
10653 ± 2 ± 2 
16.4± 3.6+g 

O.QliO.lllHa 
-20± 181^"* 


10608 ± 2l^ 
12.2 ± 1.7 ±4 
10652 ± 2 ± 2 

10.9 ± 2.6^2 
0.73±0.10lEJJ^ 
6±24t23- 


10605.1 ± 2.2+^0 
„ 4+4.5 +2.1 ■ 
ll-"* -3.9 -1.2 

10654.5 ± 2.5 +i g 
20 9+5-4+2.1 ■ 

1 o +1.0 +0-1 

-0-7 -0.5 

1 HS +44 +4 
18» _9 


10596 ±7+2 

-TO +16 +13 
^0 -10 -4 

10651 ± 4± 2 

-19+11 +8 

-, q+3-1 +0-4 

^■•^ -1.1 -0.7 
ncc +56 +12 
zoo -183 



Y{1S)7iV 
Y{2S)7tV 
Y(3S)7i*Ji' 
h^(1 P);iV 
h|^(2P)3iV 
Average 



Z,(10610) 



-10 10 
AM, MeV 



Z.(10650) 



-10 10 
Ar, MeV 



-10 10 
AM, MeV 



-10 10 
AT, MeV 



Figure 5: Comparison of ^^(10610) and ^^(10650) parameters obtained from different 
decay channels. The vertical dotted lines indicate B*B and B*B thresholds. 



1.5MeV/c^, r[Z;,(10650)] = 14.4 ± 3.2 MeV, where statistical and systematic errors 
are added in quadrature. 

The measured masses of these states exceed by only a few MeV/c^ the thresholds 
for the open beauty channels B*B (10604.6 MeV) and B*B* (10650.2 MeV). This 
"coincidence" can be explained by a molecular-like type of new states, i.e., their 
structure is determined by the strong interaction dynamics of the B*B and B*B 
meson pairs |12j . 

The widths of both states are similar and are of the order of 15MeV/c^. The 
Z;,(10610) production rate is similar to the ^^(10650) production rate for every decay 
channel. Their relative phase is consistent with zero for the final states with the 
T{nS) and consistent with 180 degree for the final states with the hb{mP). 

The T{5S) — )■ hb{mP)7i'^TT~ decays seem to be saturated by the ^^(10610) and 
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Zb( 10650) intermediate states; this decay mechanism is responsible for the high rate 
of the T(55') — )■ hh{mP)7i^7i^ process measured recently by the Belle Collaboration. 

Analysis of angular distributions for charged pions [6] favors the = 1+ spin- 
parity assignment for both Zb(10610) and Zf,(10650). Since the T(55') has negative 
G-parity, Zf, states will have opposite G-parity due to emission of the pion. 
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